Localization-induced Griffiths phase of disordered Anderson lattices 
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We demonstrate that local density of states fluctuations in disordered Anderson lattice models 
universally lead to the emergence of non-Fermi liquid (NFL) behavior. The NFL regime appears at 
moderate disorder (W = W c ) and is characterized by power-law anomalies, e. g. C/T ~ l/T 1 - 1-0 ', 
where the exponent a varies continuously with disorder, as in other Griffiths phases. This Griffiths 
phase is not associated with the proximity to any magnetic ordering, but reflects the approach to 
a disorder-driven metal-insulator transition (MIT). Remarkably, the MIT takes place only at much 
larger disorder Wmit ~ 12W C! resulting in an extraordinarily robust NFL metallic phase. 
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The nature of the non-Fermi liquid (NFL) behavior 
observed in several heavy fermion compounds remains 
largely unresolved fl. In the cleaner systems, the prox- 
imity to a quantum critical point seems to be at the origin 
of many of the observed properties J^,|| . Exotic impurity 
models cannot be discarded [Q , though their behavior in 
concentrated systems remains ill understood ||. 

In other compounds, non-stoichiometry has prompted 
the investigation of disorder-based mechanisms. A phe- 
nomenological "Kondo disorder" model (KDM), describ- 
ing a broad distribution of Kondo temperatures Tk, has 
been successfully applied to several of these systems HQ . 
Alternatively, the formation of large clusters of magnet- 
ically ordered material within the disordered phase has 
also been proposed || . Both scenarios lead to a wide dis- 
tribution of energy scales, giving rise to similar thermo- 
dynamic anomalies and NMR response j7|. In addition, 
the predictions of the KDM prove to be consistent with a 
number of other experiments, including optical conduc- 
tivity magnetoresistance |h0| ] and dynamic neutron 
scattering measurements. 

Despite these successes, a number of basic questions 
remain unresolved, including: (1) What is the micro- 
scopic origin of the ubiquitous power law (or logarithmic) 
anomalies? (2) Can a model calculation be done, which 
can produce these power laws in a universal fashion? (3) 
Are these properties tied to the proximity to a quantum 
phase transition, and if so, which one? (4) How robust is 
the anomalous behavior with respect to the variation of 
materials parameters? 

Within our model, all these questions find clear-cut 
and physically transparent answers: (i) The anomalies 
can be ascribed to a power law distribution of Tr-'s, whose 
exponent varies continuously with disorder strength. The 
resulting NFL behavior, e. g. 7 = C/T ~ l/T^- a \ 
a < 1 sets in for relatively weak randomness, irrespective 
of the detailed model for disorder. This should be con- 
trasted with the KDM J6|J^] , where the occurence of NFL 
behavior requires fine-tuning, (ii) We find universal be- 
havior reflecting the nonlocal, many-body processes asso- 
ciated with Anderson localization effects in the presence 



of strong electron correlations, (iii) for stronger disorder, 
the NFL metallic behavior persists over a surprisingly 
large interval before a disorder-driven MIT is reached. 
This novel Griffiths phase is a manifestation of quan- 
tum critical behavior associated with the approach to a 
disorder-driven metal-insulator transition and does not 
require the proximity of any magnetically ordered phase. 

We consider a disordered infinite-U Anderson lattice 
Hamiltonian 

H = E H« + e **«) C U> + E E fifUia 
ija jo 



The infinite-U constraint at each f- 



in usual notation. 

orbital is assumed (nj = fjafja — !)■ We have stud- 
ied different types of disorder, including randomness in 
the conduction electron site energies £i, the f-electron en- 
ergies Efj, as well as the hybridization Vj. Within our 
approach, we find that most of our conclusions remain 
valid for any specific form of disorder, indicating robust 
and universal behavior. 

We treat the above Hamiltonian within the recently 
proposed statistical dynamical mean field theory [fl3|| . 
This approach reduces to the usual dynamical mean field 
theory in the limit z — > 00 (with tj% ~ t/y/z) |]^, |l4] , |l5f , 
but unlike the latter, it incorporates Anderson localiza- 
tion effects. As a result, the spectral function of the local 
bath "seen" by each impurity has strong spatial fluctua- 
tions and contains information coming from sites which 
are many lattice parameters away. Physically, the fluc- 
tuations of the conduction electron wave-functions lead 
to the distribution of Kondo temperatures, which in turn 
creates a renormalized effective disorder seen by the con- 
duction electrons. This nonlocal feedback mechanism re- 
sults in the universal form of all the relevant distribution 
functions that we find. 

The simplest model for incorporating localization ef- 
fects is obtained by focusing on a Bethe lattice of coordi- 
nation z (with nearest neighbor hopping t, used as unit 
of energy). The resulting set of self-consistent stochastic 
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equations is governed by the local actions p3|Jl 

S eB = [ dr^ fUr) {dr + E fj ) /^(r) + 

[ 'dr fdr^/^WA^T-r'J^tr'); (2) 
Jo Jo a 

AjH = JL - n znr, : ■ (3) 



Here, (u>) is the local c-electron Green's function on 
site k with the nearest neighbor site j removed. It is 
determined recursively from |13|,|l6|] 
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The local self-energy E /, (ui) is obtained from the solution 
of the effective action (g) |T^,|l6|,|l^| . In order to solve the 
impurity problems, we have used the large-N mean-field 
theory at T = @. We have solved Eqs. (§§) nu- 
merically by sampling. In implementing this procedure, 
we have carried out large-scale simulations for z = 3, 
with ensembles containing up to N s = 200 sites, and 
frequency meshes containing up to N ornega — 8, 000 fre- 
quencies. The numerical integrations needed to solve the 
impurity problems have been done by a combination of 
spline interpolations and adaptive quadrature routines. 
These careful numerics have made it possible to obtain 
Kondo temperatures spanning fifteen orders of magni- 
tude, which was crucial in order to examine the long tails 
of the relevant distribution functions. 

One of the greatest advantages of our approach is its 
ability to focus on full distribution functions, which is 
essential for characterizing any Griffiths phase. Some 
typical results are presented in Fig. 0, where we show 
the evolution of the distribution of local Kondo tempera- 
tures as a function of disorder, from which one computes 
the overall response of the lattice system (See the dis- 
cussion in the first ref. of 0). We find that (Fig. 0(a)) 
the distribution has a universal log-normal form for weak 
disorder. We have verified that such a log-normal behav- 
ior is obtained irrespective of the type and shape of the 
bare disorder distribution, as long as it is not too large. 

As the disorder is increased, the distribution P(Tr-) 
no longer retains its log-normal form. Instead, a long tail 
emerges on the Iow-Tr- side, with a power law asymptotic 
form (Fig. 0(b)) 



P{T K ) ~ T t 



(a-l) 
K 



(6) 



The exponent a varies continuously with disorder, as seen 
on a plot of log(P(log(T K ))) in Fig. || 
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FIG. 1. (a) Distribution of log(TA-) for various values of 
disorder strength W/t. (b) Distribution of Tk showing the 
emergence of NFL behavior. Here, gj's are distributed uni- 
formly with width W and we have used z — 3, Ef — — 1, 
V — 0.5 and /i = —0.5, in units of t. 

Note that the value a = 1, (Figs. 0(b) and |) 
with P(Tk) ~ const., corresponds to the condition for 
Marginal Fermi Liquid behavior observed in some Kondo 
alloys j^,0] , with logarithmically divergent magnetic sus- 
ceptibility x{T) and specific heat coefficient 7. This di- 
vergent behavior becomes more singular as the disorder 
is increased past this marginal case. For example, if we 
use the simple Wilson interpolation formula for x(T) |nj 
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We thus have power law divergences with exponents 
which vary continuously with the disorder strength. If 
we take t ~ lO^K, this should be observed below a few 
tens of Kelvin. Such generic behavior has been fitted to 
some NFL compounds |l^]. Besides, x(0) will diverge at 
a critical disorder strength (W c w 1 in Fig. 0) 



X(0) 



P(Tk) 



T, 



dT, 



K 



K 



1 



(a-2) 



K 



dT, 



K 



1 W c — W 



(8) 



with a similar result for 7. Note, however, that other 
higher order correlation functions, such as the non-linear 
susceptibility X3(0); which probes higher negative mo- 
ments of the distribution (xa(0) ~ 1/TS-), will begin to 
diverge at different critical values of disorder, 



X3(0) 
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where W c z ~ 0.66 for the parameters of Fig. |^. 

This general behavior is characteristic of Griffiths 
phases pf| and should not be confused with a true phase 
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FIG. 2. Power law asymptotics of P(log(Tif)) as the disor- 
der increases. The linear behavior for small log(Tif) implies 
a power law dependence of P(Tk). Inset: the exponent a of 
Eq. (|(|). Same parameters as in Fig. |l| 

transition. The system should be viewed as a disor- 
dered metal with embedded clusters of Anderson insula- 
tors. It is precisely these poor conducting regions, with 
depleted densities of states, which give rise to imperfectly 
quenched spins and the corresponding singular thermo- 
dynamic properties. 

We should also stress that the main mechanism that 
dominates the Griffiths phase is qualitatively different 
from the one in the KDM. There, Tk fluctuations were 
simply caused by the distribution of local parameters (Vj, 
Efj) and the conduction electron DOS does not fluctuate. 
By contrast, in the present treatment, fluctuations in the 
latter are dominant. To illustrate this, all the results we 
present are obtained for a model with conduction band 
disorder only, although similar results follow for any form 
of disorder. We stress that, in a KDM treatment of this 
case, Tk fluctuations are severely limited. Here, however, 
Tk fluctuations are enhanced by the fluctuations in the 
local conduction DOS, reflecting the localization effects 
and the approach to a disorder-driven MIT. 

To confirm this picture, we examine the localiza- 
tion properties of the conduction electrons. We focus 
on the typical DOS pt vp = exp{< lnpj >}; pj — 
(l/Tr)ImG C j(uj = 0), as shown in Fig. J3. This quantity 
vanishes at any disorder-driven MIT p3| , and thus can 
serve as an order parameter for localization. Remarkably, 
we find a strong decrease of this quantity upon entering 
the Griffiths phase (W/t m 1), reflecting the strongly 
enhanced conduction electron scattering due to Kondo 
disorder. Yet, the actual localization transition, where 
the typical DOS vanishes, occurs only at much larger 
disorder (W/t rs 12). This results in a very extended 
NFL metallic region, where the thermodynamics is sin- 
gular, and the conduction electrons are almost, but not 



FIG. 3. Localization properties of the conduction electrons 
as monitored by the typical DOS as a function of disorder. 
The same quantity in the absence of Kondo spins is shown for 
comparison. The vertical dashed line indicates the boundary 
of the NFL phase. Same parameters as in Fig. [[]. 

completely localized. 

This dramatic effect has a simple physical origin. Con- 
sider the distribution of the effective scattering poten- 
tials of the conduction electrons $>j(u> — 0) (see Eq. (5)) 
introduced by the f-sites. Note that m $j(uJ = 0) = 
— ZjV 2 /£fj, where Zj is the quasiparticle weight and 
Etj the (renormalized) energy of the Kondo resonance 
at site j. For sufficient disorder, the Kondo resonances 
are randomly shifted up or down in energy, giving rise 
to $j's that can be random in magnitude but also in 
sign. The resulting distribution for the inverse quantity 
$~ 1 is shown in Fig. [| and is found to broaden with 
disorder. For W/t w 1.5, a finite density of QJ 1 = 0+ 
(i. e. <f>j = +oo) sites emerges. This is crucial, since 
the corresponding f-sites act as unitary scatterers (US's), 
characterized by a maximally allowed scattering phase 
shift (8 = tt/2) for the conduction electrons. If all the 
f-sites were US's, the system would be a Kondo insula- 
tor. The presence of a finite fraction of US's should be 
viewed as the emergence of droplets of a Kondo insula- 
tor within the heavy metal. Interestingly, at stronger 
disorder (W/t > 4) the distribution of $~ 1 (0) contin- 
ues to broaden, leading to a decrease in the number of 
US's. This is illustrated by plotting P^J 1 = 0) in the 
inset of Fig ^. In this regime, while the bare disorder 
increases, the effective disorder produced by the f-sites is 
reduced, stabilizing the almost localized metallic phase. 
This mechanism may be at the origin of the puzzling 
behavior of materials such as SmBg |2lJ], where the low 
temperature resistivity remains anomalously large yet fi- 
nite over a broad range of parameters. 

Finally, we note that a similar NFL phase was iden- 
tified in a study of the Mott- Anderson transition in the 
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FIG. 4. Distribution function for the quantity $~ 1 (0) as 
a function of disorder. The inset shows the behavior of 
the concentration of the strong scattering Kondo centers 
( ( J>^ 1 (0) = 0), which reaches a maximum just after the NFL 
phase is entered. Same parameters as in Fig. ^| 



disordered Hubbard model |K|. We have re-examined 
this system, and concluded that the NFL behavior should 
be attributed to a related Griffiths phase rather than a 
separate thermodynamic phase of the system. Despite 
the similarities, several features prove dramatically dif- 
ferent. For Hubbard models, the emergence of NFL be- 
havior does not have a dramatic effect on the conduction 
electrons and no US's emerge. This observation may ex- 
plain the strong correlation between thermodynamic and 
transport anomalies in Kondo alloys, but not in doped 
semiconductors. In the latter materials, the thermody- 
namics is still singular close to the MIT, while transport 
remains more conventional |2^ ]. 

It would be of particular interest if it could be tested 
experimentally whether these localization effects are re- 
sponsible for the observed NFL behavior of disordered 
heavy fermion systems. A scanning tunneling microscopy 
study might be able to detect the predicted insulator 
droplets. In order to distinguish this from the magnetic 
Griffiths phase scenario [||, a systematic study of sys- 
tems with comparable amounts of disorder but different 
magnetic character would be useful. Besides, since the 
present theory relies very little on intersite magnetic cor- 
relations, a determination of the typical size of the rele- 
vant magnetic clusters could also serve as a test. 

In summary, we have investigated and solved a micro- 
scopic model for disordered Anderson lattices that dis- 
plays an unprecedented sensitivity to disorder, leading 
to localization-induced non-Fermi liquid behavior. Our 
results demonstrate that a well defined Griffiths phase 
can exist, which is not restricted to the vicinity of any 
magnetic ordering and yet seems to be consistent with 
most puzzling features of disordered heavy fermion sys- 



tems and Kondo alloys. 
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